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A spatially extended classical system with metastable states subject to weak spatiotemporal 
noise can exhibit a transition in its activation behavior when one or more external parameters are 
varied. Depending on the potential, the transition can be first or second-order, but there exists 
no systematic theory of the relation between the order of the transition and the shape of the 
potential barrier. In this paper, we address that question in detail for a general class of systems 
whose order parameter is describable by a classical field that can vary both in space and time, and 
whose zero-noise dynamics are governed by a smooth polynomial potential. We show that a quartic 
potential barrier can only have second-order transitions, confirming an earlier conjecture [l[ . We then 
derive, through a combination of analytical and numerical arguments, both necessary conditions and 
sufficient conditions to have a first-order vs. a second-order transition in noise-induced activation 
behavior, for a large class of systems with smooth polynomial potentials of arbitrary order. We find 
in particular that the order of the transition is especially sensitive to the potential behavior near 
the top of the barrier. 



INTRODUCTION 

When a spatially extended classical system with multi- 
ple locally stable states is perturbed by weak spatiotem- 
poral noise, a transition in its activation behavior can 
occur as one or more parameters of the system are var- 
ied [3, 0, S] • In the simplest one-dimensional systems this 
parameter is simply the length of the interval on which 
the system is defined, but for more complicated systems 
other parameters come into play. For example, a transi- 
tion from Arrhenius to non-Arrhenius behavior in ther- 
mally activated magnetization reversal in thin annular 
nanomagnets can occur either as ring size increases or as 
the externally applied magnetic field decreases [1, [j| . An- 
other example is a crossover from uniform to instanton- 
like decay of metastable metal nanowires [6| as either the 
length of the nanowire or the stress applied to it is varied. 
A similar crossover, from thermal activation to quantum 
tunneling, occurs in various systems as t emp erature is 
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These two cases are formally related through a map- 
ping that identifies interval length (and/or magnetic field, 
stress, if appropriate) in the classical field case to tem- 
perature in the quantum case; in particular, increasing 
interval length in the former corresponds to the lowering 
of temperature in the latter, with thermal activation in 
a classical field theory of infinite domain size mapping 
to zero-temperature tunneling in quantum field theories. 
These transitions are fundamentally different from the 
more usual sort, in which a change in order parameter 
(i.e., expectation value of the field) results from vary- 
ing a control parameter, such as coupling strengths in 
the Hamiltonian or noise amplitude. (For an extensive 
discussion of this more conventional kind of transition 



within a noisy field-theoretical framework, see, for exam- 
ple, [H}.) The extent to which the more unconventional 
type of transition under discussion here can be compared 
to a true second-order phase transition, and where the 
analogy breaks down, was discussed extensively in [l9| . 
(For recent work on discretized versions of these and sim- 
ilar models, see 2ol 2lj|.) 

Chudnovsky [3| first noted that the 
classical^>quantum transition in an extended sys- 
tem with a single degree of freedom can be either first- 
or second-order, depending on the potential. This 
observation has important applications. First-order 
transitions in the classical-quantum escape rate have 
been considered in anisotropic bistable large-spin mod- 
els [iil and biaxial spin systems subject to longitudinal 
fields [23], and have been more generally considered 
in the decay of metastable states in quantum field 
theories A general discussion can be found 

in El. 



Similarly, a first-order transition can also occur in 
classical transitions between two thermally activated 
regimes, as was recently found by the authors Q in the 
decay of nanowires due to thermal fluctuations (cf. Fig. 2 
of Ref. Q). Despite its potential importance, relatively 
little systematic work has been done to identify the gen- 
eral conditions under which one or the other kind of 
transition occurs. It was conjectured [l| that smooth 
polynomial potentials with terms no higher than quartic 
display only second-order transitions. It is also known 
that hi ghe r-order terms can lead to a first-order tran- 
sition [24j. At the present time, however, there is no 
systematic theory of how the order of the transition de- 
pends on potential characteristics. The purpose of this 
paper is to address that problem. 
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MODEL 

We will consider a general class of models of extended 
systems describable by a classical field <fi(z,t) defined on 
the spatial interval [-L/2, L/2], subject to a potential 
V(<p) and perturbed by spatiotemporal white noise [26j |. 
Time evolution is governed by the stochastic Ginzburg- 
Landau equation 



a t (t> = d z A - a+vw + V2T£(z, t) 



(i) 



where all dimensional quantities have been scaled out. 
The first term on the RHS arises from a field 'stiff- 
ness', i.e., an energy penalty for spatial variations of 
the field. The noise t) satisfies (£(2i,ii)£(22,*2)) = 
S(zi — Z2)5(ti — £2), and its magnitude T is small com- 
pared to all other energy scales in the problem (formally, 
our analysis will be asymptotically valid in the T — > 
limit) . 

The zero-noise dynamics of ([1]) can be written as the 
variation of an action TL with the field d>: 



d t( j> = -mm 



with 



HW\ = 



L/2 



L/2 



V{4>) 



(2) 



(3) 



Stationary solutions of @ describe stable, metastable, 
and transition (i.e., saddle) states of the system. 

In the weak-noise (T — > 0) limit, the classical activa- 
tion rate for a transition out of a (meta)stable well is 



T - T exp(-AW/T) 



(4) 



where the activation barrier AW is the action difference 
between the transition state and the initial (meta)stable 
state. The rate prefactor T is determined by fluctuations 
about the most probable escape path. When the top 
of the barrier is locally quadratic, the prefactor Tq is 
independent of temperature. In such circumstances the 
escape rate Q is said to be of the Arrhenius-van't Hoff 
(or often simply Arrhenius) form. Here we will mostly 
be concerned with the behavior of the activation barrier 
AW. 



TRANSITION IN ACTIVATION BEHAVIOR 

We briefly summarize here the derivation of a second- 
order transition in the noise-induced barrier crossing de- 
scribed by |T]) between wells in the simple bistable sym- 
metric quartic potential 



V s {4>) = (4? - l) 2 , 

with Neumann boundary conditions dz£ 
Qz4\l/2 — 0- The discussion follows that of |l| 
we refer the reader for details. 



(5) 
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FIG. 1: (a) The activation energy AW and (b) the rate 
prefactor To as functions of the interval length L, for the po- 
tential given by Eq. (JSJ) with Neumann boundary conditions. 
The dashed line indicates the critical interval length L c = n/2 
at which the saddle state bifurcation takes place, showing the 
power-law divergence of the prefactor. The transition state 
4>t(z) for L — 5 (corresponding to m = 0.986) described by 
Eq. (JSJ is displayed in the inset (only one of the symmetric 
pair is shown.) Note that quantities in all figures are expressed 
in dimensionless units (see text.) 



Because of the symmetry of the potential (which is 
not necessary for the transition to occur [H), the change 
in activation behavior arises from a bifurcation of the 
transition state. Below a critical length L c the transition 
state <f>t is constant, while above L c it becomes a pair of 
degenerate, spatially varying instanton configurations 0]: 



0. 



±\/i 2 r 2 -sn( 

V l+m V 



L <L C 
L> L r 



(6) 



where sn(- | m) is the Jacobi elliptic sn function with pa- 
rameter < m < 1. Its quarter-period is given by K(m), 
the complete elliptic integral of the first kind 27[ , which 
is a monotonically increasing function of m. As m — > + , 
K(m) decreases to n/2, and sn(- | to) — * sin(-). In this 
limit the saddle state smoothly degenerates to the <j> t = 
configuration. As m — * I~, the quarter-period increases 
to infinity (with a logarithmic divergence), and sn(- | 
to) — > tanh(-), the (nonperiodic) single-kink sigmoidal 
function. The Langer-Callan-Coleman [H, Hi,B3, 
'bounce' solution is thereby recovered as L — > 00. 

The value of m in © is determined by the interval 
length L and the Neumann boundary conditions, which 
require that 



L = \Jm + I K(to) 



(7) 



The critical length is determined by when m = 0; that 
is, L c = n/2. As previously noted, m — ► 1 corresponds to 
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FIG. 2: Activation barrier AW (black lines) for various po- 
tentials exhibiting first-order transitions, marked by open cir- 
cles. The line styles correspond to the following potentials: 
(solid) V(<p) = 1 - #73 - 4 + 4073, (dashed) V{4>) = 
1-0 2 — ^ 3 +0 5 , (dot-dashed) V{4>) = 1-2073-3072+7076, 
and (dotted) 1/(0) = 1 - 130712 - 3 + 074 + 5076, with 
the curves shifted horizontally to improve readability. The 
(second-order) transition between uniform and instanton-type 
escape is marked by an open square for each curve. Activation 
barriers for higher energy states are shown with gray lines. 
Degenerate instanton-like saddle states corresponding to the 
first-order transition for the dotted line are shown in the in- 
set. Both states have the same length and energy. Note that 
the potential described by the rightmost curve (dotted line) 
exhibits a second-order followed by a first-order transition. 



L — ► oo, and the activation energy smoothly approaches 
the asymptotic value of AWoo = 4v2/3. The transition 
state for an intermediate value of m, corresponding to 
L = 5, is shown in the inset of Fig. [T] 

The activation energy AW can be computed in closed 
form for all L > L c (below L c , it is simply AW = L): 



AW 



[1 + m)E(m) - (1 - m)(3m + 5)K(m) 



3(1 



\3/2 



(8) 



with E(m) the complete elliptic integral of the second 
kind [271 ]. The activation energy as a function of L is 
shown in Fig. QJa). Note that the curve of AW vs. L 
and its first derivative are both continuous at L c ; the 
second derivative, however, is discontinuous, as might be 
expected of a second-order-like phase transition. 

A more profound manifestation of critical behavior at 
L c is exhibited by the rate prefactor Tq, which (in the 
asymptotic limit T — > 0) diverges at L c , as shown in 
Fig. [Hb). This is striking, but requires interpretation. 
Because it is not relevant to the present discussion, we 
refer the interested reader to 

This behavior is generic for a whole class of poten- 
tials, as described below, but first-order transitions have 
also been observed jfjj, leading to a continuous but non- 
differentiable activation barrier AW at the transition 



points. This is illustrated for various potentials in Fig. [2] 
The second-order transition is still present, but is usually 
not physically observable, as it happens for higher energy 
transition states (note, however, the exception in Fig. 2, 
dotted line). 



ORDER OF THE TRANSITION 

Eqs. (|2|) and (J3j) lead to the (typically nonlinear) dif- 
ferential equation for stationary states 



0" = F'(0). 



(9) 



If we map the field to position and the coordinate z 
to time, the solutions to this equation are equivalent to 
the trajectories of a classical particle moving in the in- 
verted potential — V(4>) [32]. The bounce, or 'instanton', 
transition state (cf. Eq. ^ that determines the activation 
behavior when L > L c ) then corresponds to a half-period 
of such a periodic classical trajectory. 

These classical trajectories have a corresponding "en- 
ergy" E = -E(L) given by (cf. Eq. (6) of Q) 



E{L) = \{<t>'?-V{cl>). 



(10) 



This energy of a classical instanton trajectory should not 
be confused with the activation barrier AW given by the 
action difference between the transition and metastable 
states. E corresponds to the energy of a classical parti- 
cle undergoing periodic motion in the inverted potential 
— V(<p). It is determined either by the temperature in 




E 



FIG. 3: The energy of the classical trajectory determining 
activation behavior for (solid line) the symmetric quartic po- 
tential of Eq. (0), (dashed line) v\<j>) = 1 - 50 72 + 3 + 4 /2, 

- 20 2 + 072 + 072, and (dotted 



(dot-dashed line) V(<j>) = 1 ■ 

line) V(0) = 1 - 7074 + 073 + 50 6 /12. The corresponding 
potentials V(<j>) are plotted in the inset. As discussed in the 
text, the monotonic decrease of energy with interval length 
signifies a second-order phase transition at L c . 
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FIG. 4: The energy of the classical trajectory determining ac- 
tivation behavior for various potentials exhibiting first-order 
transitions, as evidenced by the non-monotonic behavior of 
L(E). The line styles correspond to those of Fig. [2] with the 
corresponding potentials V(4>) plotted in the inset. 



the thermally assisted tunneling problem, or by the in- 
terval length in a stochastic classical Ginzburg-Landau 
field theory. 

We illustrate the behavior of the energy for the simple 
case of the symmetric quartic potential given by Eq. ([5]). 
For L < L c , the transition state 0j = [cf. Eq. ©], and 
the energy E = — E = 1, independent of length. For 
L > L c , the energy E monotonically decreases to zero as 
interval length increases. As a function of to (related to 
L through Eq. 0), it is given by 



E(m) = 1 - 



4m 



(1 + m) 2 



(11) 



The length as a function of energy is shown in Fig. [3] 
for various potentials, plotted in the figure inset, with 
the solid line corresponding the Eq. ([5]) . It is instructive 
to compare this to the behavior of the activation energy 
AW as a function of length (Fig. [1]). 

As first noted by Chudnovsky [14|, the order of the 
transition is related to behavior of the period L of the 
instanton trajectory vs. its energy E: a monotonic de- 
crease, as in Fig. [3J signifies a second-order transition 
while a non-monotonic decrease corresponds to a first- 
order transition, as shown in Fig. [3] for a variety of po- 
tentials. 

Our interest here is in determining the relation between 
the potential properties and the order of the transition. 
In order to facilitate comparison between different po- 
tential barriers, we rescale 7(0) to unit barrier height 
and width, the latter being defined as the distance be- 
tween the maximum and minimum of 7. Specifically, we 
rescale V(<f>) so that 7(0) = 1, 7'(0) = 0, 7"(0) < 
and 7(1) = 0, 7'(1) = 0, 7"(1) > 0. We will generally 
take the state = 1 to be metastable; it can decay to- 



ward negative values of if there is a e < such that 
7(0 e ) < 0. 

We begin with some general considerations, valid for 
all potential barriers. All L(E) curves in both Figs. [3J 
and S] diverge as E — > + . This generic behavior is easily 
understood from the classical-particle analogy, as a par- 
ticle with a periodic orbit will have an arbitrarily long 
period if its energy is close to, but lower than, a local 
maximum of the potential —7(0). 

Similarly, a barrier containing a local, secondary maxi- 
mum/minimum pair (see, e.g., the dot-dashed line in the 
inset of Fig. [4} will have a divergence of L(E) at the en- 
ergy corresponding to the secondary minimum. Thus, a 
potential barrier that includes a local metastable state 
leads to a first-order transition of the activation behav- 
ior. That transition, in this case, is between a two-step 
escape through the local minimum and a direct escape, 
with both paths having instanton- like transition states. 

As L{E) initially decreases for increasing E when E is 
small, a sufficient condition for having a first-order tran- 
sition is that dL/dE > for E — ► 1 _ . The corresponding 
condition on the potential 7 can be derived analytically 
using perturbation theory around the barrier maximum. 

More general conclusions require a numerical determi- 
nation of L{E). To simplify the analysis, but still keep 
the conclusions general, we consider smooth potentials of 
the form 



7(0) = 1 - a 2 2 + a m m + i 



(12) 



with 2 < m < n. Rescaling to unit barrier height and 
width leaves a single free parameter —nj (n — m) < a m < 
2/ (to - 2), with 



a 2 



and 



n + (n — m)a r , 



(to — 2)a r , 
n-2 



(13) 



(14) 



Our approach is to solve numerically the nonlinear dif- 
ferential equation corresponding to a particle in an in- 
verted potential —7(0) with initial condition 0(— L/2) = 
0o, with < 0o < 1, and 0'(— L/2) = 0, and with L the 
minimal length satisfying Neumann boundary conditions 
at z = L/2. i?(0o) follows from Eq. (|10p . thus providing 
the function L{E) in parametric form. 

Our findings are summarized graphically in Fig. 
where — 2 /2 + 0™/n (solid lines) and a<p m /m (dashed 
lines) are plotted separately. Potentials having first-order 
transitions are shown on the left-hand side of the dotted 
line, and those with second-order transitions are on the 
right-hand side. All potentials on the left-hand-side have 
a negative m term, which appears to be a necessary, al- 
though not a sufficient, condition for the existence of a 
first-order transition. 
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FIG. 5: (Color online) Var- 
ious types of potentials of 
the form given by Eq. (|12[) . 
with the second term a<f> m /m 
plotted as a (blue) dashed 
line, and the other terms, 
-0 2 /2 + (f) n /n, plotted as a 
(red) solid line. Potentials on 
the left-hand-side of the dot- 
ted line have first-order tran- 
sitions, while those on the 
right-hand-side have second- 
order transitions. 



More specifically, for a potential of the form (fT2|) . the 
transition is unique and second-order if either of the 
following conditions is fullfilled: 

• a m Ss 0; 

• a m < 0, with m = 2k + 1 and n = 21, k and I being 
positive integers such that 1 < k < I. 

The second condition is actually a subset of the first 
for the symmetric potential V(—<p). This confirms in par- 
ticular that polynomial potentials with at most quartic 
terms exhibit only second-order transitions [H. 

On the other hand, the transition is first-order if 
a m < and 

• m = 2k, with k > 1, 
or 

• m = 2k + 1 and n = 21 + 1, with 1 ^ k < I. 

If the third term on the RHS of Eq. (fT2|) is an even 
power of <j), then the transition is first-order if and only 
if a m < 0. If it is an odd power of <j>, the transition is 
first-order if the third term is an odd power of <fi as well, 
but with an opposite sign, so that there is a competition 
between the two terms. The mechanism leading to a first 
order transition thus seems to be different depending on 
the parity of the third term in V, as discussed below. 

POTENTIALS WITH A SECOND-ORDER 
TRANSITION 

Potentials given by Eq. (JT2J) have a second-order tran- 
sition in the following cases: 



• m = 2k, for k > 1, and a m ^ 0; 

• m = 2k + 1, for k ^ 1, and at least one of the 
conditions n — 21 or a m 0. 

These potential barriers all look very similar, and have 
the same generic behavior, shown in the inset of Fig. [3] 
for various combinations of m, n, and a m : The function 
L(E) decreases monotonically for E > 0, and the escape 
energy has a linear part corresponding to a homogeneous 
transition for L < L c , with a second-order transition to 
instanton-like escape at L c = 7r/y2a2. 

The transition length is set by the — <fi 2 term in the 
potential: a potential V{4>) = 1 — gives an equa- 

tion for the transition state that is equivalent to a har- 
monic oscillator in one dimension, and therefore has a 
transition state whose length L = ir/ \J2ot2 and energy 
E u = 1 are independent of the initial value <po ■ The term 
proportional to <f> n modifies that behavior by increasing 
the length of the transition state for a given initial value 
00 compared to the inverted quadratic potential. In the 
analogy to a particle in a one-dimensional potential —V, 
it decreases the slope of the potential for larger <j), propor- 
tional to the force felt by the particle, thus increasing the 
period of the orbit, until one approaches the maximum 
of —V, where the orbital period diverges. 

Considering the case a m — 0, one can study numeri- 
cally the ratio of the escape barrier A Woo for large inter- 
val length L to the energy at the transition point AW(L C ) 
as a function of n (see Fig. |6|). This ratio approaches 1 
as a power law with exponent 3/2 for n 3> 1, with some 
even-odd oscillations. As n becomes large, the transition 
therefore looks more and more like a first-order transi- 
tion, especially for odd n. Additionally, the curve L(E) 
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FIG. 6: Ratio of the saturation energy barrier A Woo to the 
energy barrier at the transition point AW(L C ) as a function 
of the power n of the last term in the potential V{4>) — 1 — 
^2 2 + ^2 0™! with circles (triangles) marking odd (even) 
values of n. The lines are a guide-to-the-eye showing power- 
laws with exponenet 3/2. 



case, thus decreasing the oscillation period. As E — > 1~, 
the quadratic term dominates and limits the period to its 
value at L c . As the energy E decreases the last term in 
Eq. (fT2")l inverts the trend and increases the period again, 
thus creating a minimum in L(E). 

This result can be generalized to the potential 



v(<t>) 



i 



h2 E 

k>0 



a k <j) k 



(16) 



with m > and ao ^ 0, for energies E(e) = 1 — uj 2 e 2 , 
e -C 1, using perturbation theory. Using Eq. (fTU|) . one 
can derive an expression for the half-period 



L(E) 



W _e) ^2[V{4>)-E] 



(17) 



where <X±e) is the smallest on the right/left-hand- 
side of the maximum such that V(</>(±e)) = E(e). Using 
a series expansion in e for </>(±e), and expanding Eq. |(T7|) 
to order e 2m , one obtains 



2,2^ 



remains almost flat for an increasing range of E < 1 as n 
becomes larger. This implies that, at the critical length 
L c , there is a continuum of transition states with quasi- 
degenerate energies available for the escape. 



POTENTIAL WITH A FIRST-ORDER 
TRANSITION: EVEN m = 2k 

Negative terms of power m larger than 2 in the Taylor 
expansion of V((j>) around its maximum are responsible 
for the change of order of the transition, but the mech- 
anism creating a first-order transition depends on the 
parity of m. 

In the case of an even m — 2k > 2, with a m — — |o!2fe|, 
the function L(E) has a local minimum at some energy 
< E < 1. The solid line in Fig. 0] shows a typical 
behavior for that type of potential. The final increase 
of L(E) as E — ► 1, corresponding to an initial decrease 
of L(4>q) for <fio <C 1, is driven by the middle term in V, 
— | c»2k 1 4> 2k '■ on its own, such a term creates a divergence of 
the curve L(E) for E — > 1, or </>o — + 0, with L thus being 
an increasing function of E. One can derive an analytic 
expression for L(E) for a potential V(tf>) = 1 — a4> m /m 
by using Eq. (fTU|) to express dz in terms of d<p and V(^). 
Integrating that expression yields 



L(E= 1-u/e 



L(E) = 



2irT{l/m)a- 1 / m 
;T[(m + 2) /2m] 



(1-E)^ 



(15) 



The presence of the quadratic term removes the diver- 
gence, but the increase for L(E) as E — > 1 remains. In 
terms of the classical particle, the — (f> 2k term increases 
the initial force on the particle compared to the quadratic 



fc=0 



1 + (-l) T " +fc ^k A , c m +k 

2 



JlB 2m e 2m + 0(e) 



2m+l 



(18) 



where the numbers Ak > 0, k — 1, . . . m and B2m > 0. 
From Eq. (JTHJ) it is clear that dL/dE = -l/(2u> 2 e) • 
dL/de > for £7 > 1 if the lowest even-power 2n > 2 
term in the potential V(cj>) has a negative coefficient, re- 
gardless of the presence of odd-power terms, which do not 
contribute to the slope of L(E) to that order, provided 
n < m. 

This provides a sufficient, though not necessary, condi- 
tion for the presence of a first-order transition: the tran- 
sition is first-order if the first even power of cj> (excluding 
the quadratic term) in the Taylor expansion of the po- 
tential around its maximum has a negative coefficient, 
provided its exponent is at most 2m, with m defined in 
Eq. CBJ. 



POTENTIAL WITH A FIRST-ORDER 
TRANSITION: ODD m = 2k + 1 

Perturbation theory yields no information about the 
influence of odd-power terms in the potential as they do 
not contribute (up to order 0(e) 2m ) to the slope of L(E) 
for E — * 1~. However even a negative slope does not 
imply a second-order transition, as illustrated in Fig. [3] 
(dotted and dashed lines). In order to study the effect of 
odd-power terms, we return to numerics, with a potential 



V(4>) = 1 - a 2 r - \a k \(j) Ik+l + 



k 2t+l 



(19) 
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with 1 < k < I. 

In this case, the function L(E) decreases for E — ► 1~, 
but has both a maximum and a minimum for lower ener- 
gies. For |offc| > a cr it, the maximum of L((po) turns into 
a divergence. This is related to the formation of a sec- 
ondary maximum in the potential V, as is the case for the 
dot-dashed line in Fig. [5J where ctk is slightly larger than 
the critical value for that particular potential. The max- 
imum of L(E) for ctk < a cr it comes from the 'formation' 
of that secondary maximum with increasing at- 

Note that the transition from homogeneous to 
instanton-like escape (marked by squares in Fig. [51 hap- 
pens above the lowest activation energy curve (black line) 
for potentials given by Eq. (fT2|) . and is therefore not ob- 
servable. However the addition of a positive quartic term 
to the potential, as is the case for the dot-dashed line 
of Figs. [2] and 0J can move that second-order transition 
'below' the first-order one. In that case there would be 
two transitions: a second-order transition from a homo- 
geneous escape to escape through an instanton, followed 
by a first-order transition between two different escape 
routes with instanton transition states (shown in the in- 
set of Fig. . 



DISCUSSION 

We have presented a comprehensive study of the de- 
pendence of the order of the barrier crossing transition on 
potential characteristics, for classical extended systems 
subject to weak external spatiotemporal noise. Using 
a combination of analytical and numerical methods, we 
confirmed an earlier conjecture of one of the authors [lj 
that smooth potentials whose highest term is quartic 
have second-order transitions. We then considered a wide 
class of polynomial potentials of arbitrary order, and de- 
termined the potential characteristics that led to cither 



TABLE I: Summary of results: The first column lists the 
transition order (I or II) for the potential V(cj>) = —(f> 2 + F((j)), 
with F(4>) given is the second column. The third column lists 
conditions that F(<p) needs to satisfy. 



F{4>) 



Conditions 



II 
II 



j.2fc+l . ,21 
OLk<t> + 4> 



a m > 0, 2 < m < n 
a k < 

2 < 2k + 1 < 21 



oi k <p + <p 



,2k+l . ,21+1 

a k (j> — (— d) 



I 7n0 2 " + ^ Qfct 



2(m+fe) + l 



+ ■ 



a k < 0, 2 < 2k < n 
Q fe < 0, 1 < k < I 
m > 1, n > 1, n < 2m 
aso / 0, 7„ < 



first- or second-order transitions. These results are sum- 
marized in Fig. [5] and Table |TJ In particular, we found 
that the potential characteristics at the top of the bar- 
rier play a central role in determining the order of the 
transition. 

The order of the transition can play a crucial role in 
the understanding of systems near the transition point 
and in the design of new experiments (and possibly de- 
vices). For example, in [331] . a transition from ohmic to 
non-ohmic behavior was observed as the length of the 
gold nanowires was changed. This was explained 34 1 
in terms of the transition predicted in [f| for monova- 
lent metallic nanowires as wire length changes. So there 
already exists experimental evidence for a transition. A 
more systematic experimental investigation of this ohmic 
to non-ohmic transition is highly desirable in order to ex- 
amine details of the behavior near the transition point. 
Meanwhile, the analysis given here allows predictions as 
to which wires (characterized by their radius, or more 
easily measurable, long-wire conductance) will undergo 
one or the other type of transition. These can be mea- 
sured as a difference in behavior — sharp discontinu- 
ity vs. smooth crossover — of an effective temperature 
T e f f (T) characterizing the barrier crossing rate out of a 
given metastable state (22}. Other applications can be 
found in [13. 

Although our studies are done in the context of classi- 
cal transitions between metastable states, they should be 
generally applicable to a broad set of problems, including 
the classical<->quantum crossover or transitions between 
regimes of thermally-assisted quantum tunneling, follow- 
ing the mapping described in [191 ]. 

This work was supported by NSF Grant Nos. 0312028 
(CAS) and PHY-0651077 (DLS). 
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